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Abstract 

Double porosity models for the liquid filtration in a naturally 
fractured reservoir is derived from the homogenization theory. The 
governing equations on the microscopic level consist of the station- 
ary Stokes system for an incompressible viscous fluid, occupying a 
crack-pore space (liquid domain), and stationary Lame equations for 
an incompressible elastic solid skeleton, coupled with corresponding 
boundary conditions on the common boundary "solid skeleton-liquid 
domain" . We suppose that the liquid domain is a union of two inde- 
pendent systems of cracks (fissures) and pores, and that the dimen- 
sionless size 5 of pores depends on the dimensionless size e of cracks: 
6 = e"' with r > 1. The rigorous justification is fulfilled for homog- 
enization procedure as the dimensionless size of the cracks tends to 
zero, while the solid body is geometrically periodic. As the result we 
derive the well-known Biot - Terzaghi system of liquid filtration in 
poroelastic media, which consists of the usual Darcy law for the liquid 
in cracks coupled with anisotropic Lame's equation for the common 
displacements in the solid skeleton and in the liquid in pores and a 
continuity equation for the velocity of a mixture. The proofs are based 
on the method of reiterated homogenization, suggested by G. Allaire 
and M. Briane. As a consequence of the main result we derive the 
double porosity model for the filtration of the incompressible liquid in 
an absolutely rigid body. 
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Introduction 



The liquid motion in a naturally fractured reservoir is described by different 
mathematical models. These models take into account a geometry of a space, 
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occupied by the liquid (liquid domain), and physical properties of the liquid 
and the solid skeleton. Among different models the simplest one is Darcy 
equations 

v = -kVq + F, \/-v = 0, (0.1) 

for the macroscopic velocity v and the pressure q of the liquid, when the solid 
skeleton is supposed to be an absolutely rigid body and the liquid domain 
is a pore space. For more complicate geometry, when the liquid domain is 
a union of system of pores and cracks, there are different type of models 
(see, for example, Ref. Ref. HI], Ref. [20], Ref. [22]). Note, that pores 
differ from cracks by its characteristic size: if Ip is a characteristic size of 
pores and Ic is a characteristic size of cracks, then Ip -C Ic- The well-known 
double-porosity model, suggested by G. I. Barenblatt, lu. R Zheltov and I. 
N. Kochina [1], describes two- velocity continuum where macroscopic velocity 
Vp and pressure in pores and macroscopic velocity and pressure Qc in 
cracks satisfy two different Darcy laws 

Vp = -kpVqp + F, v, = ~k,Vqc + F, (0.2) 

and two continuity equations 

V ■Vp = J, V ■v^ = - J. (0.3) 

The model is completed by postulating that the overflow J from pores to 
cracks linearly depends on the difference (gc — qp)- 

In view of the importance of such models it is very natural to rigorously 
derive the governing equations for each model, starting with detailed mi- 
crostructure of the liquid domain and the linearized equations of fluid and 
solid dynamics on the microscopic level. In their fundamental paper R. Bur- 
ridge and J. Keller [8] have used this scheme to justify a well - known in 
contemporary acoustics and filtration phenomenological model of poroelas- 
ticity, suggested by M. Biot [5J. As a model of the porous medium on the 
microscopic level authors have considered the mathematical model, consisting 
of Stokes equations describing liquid motion in pores and cracks, and Lame's 
equations, describing motion of a solid skeleton. The differential equations 
in the solid skeleton and in the liquid domain are completed by boundary 
conditions on the common boundary "liquid domain - solid skeleton" , which 
express a continuity of displacements and normal tensions. The suggested 
microscopic model is a basic one, because it follows from basic laws of con- 
tinuum mechanics ( see also E. Sanchez - Palencia [19j). After scaling there 
appears a natural small parameter 6 which is the pore characteristic size Ip 
divided by the characteristic size L of the entire porous body: 6 = Ip/L. The 
small parameter enters both into coefficients of the differential equations, and 
in the geometry of the domain in consideration. The homogenization (that 
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is a finding of all limiting regimes as (5 \ 0) of this model is a model, asymp- 
totically closed to the basic model. But even this approach is too difficult 
to be realized, and some additional simplifying assumptions are necessary. 
In terms of geometrical properties of the medium, it is most expedient to 
simplify the problem by postulating that the porous structure is periodic 
with the period 6. Under this assumption R. Burridge and J. Keller, using a 
method of two-scale asymptotic expansion, have formally justified M. Biot's 
model. For the same geometry of the pore space (let call such a model as a 
single porosity model) and for absolutely rigid solid skeleton when a liquid 
motion is described by the Stokes system, L. Tartar have rigorously justified 
the Darcy law of filtration (see Appendix in Ref. [IHj). Later a rigorous jus- 
tification of M. Biot's models, under same assumptions on the geometry of a 
pore space as in Ref. [S], has been rigorously proved in Ref. [13] - Ref. P^ . 



For more complicate geometry, when the liquid domain is a crack - pore 
space (let call such a geometry as a double porosity geometry and corre- 
sponding mathematical model as a double porosity model), some attempts to 
derive macroscopic models, asymptotically closed to some phenomenological 
models on the microscopic level have been made by T. Arbogast et al [3j, 
A. Bourgeat et al [7] and Z. Chen [9]. Because the last two papers repeat 
ideas of the first one, let us briefly discuss the main idea in Ref. [3]. As a 
basic model on the microscopic level, the authors have considered a peri- 
odic structure, consisting of "solid" blocks of the size e surrounded by the 
fluid. The solid component is assumed to be already homogenized: there 
is no pore space and the motion of the fluid in blocks is governed by usual 
Darcy equations of filtration. The motion of the fluid in crack space (the 
space between "solid" blocks) is described by some artificial system, similar 
to Darcy equations of filtration. There is no any physical base, but from 
mathematical point of view such a choice of equations of fluid dynamics in 
cracks is very clear: it is impossible to find reasonable boundary conditions 
on the common boundary "solid" block-crack space, if the fluid dynamics 
is described by the Stokes equations. But there are reasonable boundary 
conditions, if the liquid motion is described by Darcy equations of filtration. 
Therefore, the final macroscopic models in Ref. [3], Ref. [7j and in Ref. [H] 
are physically incorrect (see Ref. [T7j). 

The physically correct double porosity model for the liquid filtration in 
an absolutely rigid body has been derived by A. Meirmanov p2j. Following 
the scheme, suggested by R. Burridge and J. Keller [S], author starts with 
a liquid domain, composed by a periodic system of pores with dimensionless 
size 6 and a periodic system of cracks with dimensionless size e, where 6 = e^, 
r > 1. The liquid motion is described by the Stokes system 
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arpf^ = af,^v -Vq + PfF, — + OgV ■ v = 0, 



(0.4) 
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for dimensionless microscopic velocity v and pressure q of the liquid, where 

L 2n c^pf 

gr^ TLgpo Lg 

L is a characteristic size of the domain in consideration, r is a characteristic 
time of the process, p/ is the mean dimensionless density of the liquid, scaled 
with the mean density of water pq, g is the value of acceleration of gravity, p 
is the viscosity of fluid, c is a speed of sound in fluid, and the given function 
F{x,t) is the dimensionless vector of distributed mass forces. 

It is supposed, that all dimensionless parameters depend on the small 
parameter e and the (finite or infinite) limits exist: 

limar(e) = rn, limaJe) = pn, Imxanie) = c,, lim — ^ = ui, lim — ^ = 

e\0 £\0 e\0 ^ e\0 s\0 6^ 

The aim of any homogenization procedure of some mathematical model, de- 
pending on the small parameter e, is to find all possible limiting regimes in 
this model as e \ 0. Of course, these regimes for the model (10.41) depend on 
criteria tq and pi, which characterize different types of physical processes. 
We may roughly divide all these processes on two groups: long-time processes 
(filtration) and short-time processes (acoustics). It is well-known, that the 
characteristic time of the liquid filtration is about month, while the char- 
acteristic size of the domain is about thousand meters. Therefore, we may 
assume that for filtration tq = 0. The rest of processes we call acoustics and 
all these situations characterized by criterion tq > 0. 
Under restrictions 

Po = 0, To < oo, < C/ < oo, 

the author has shown that the homogenization procedure for the liquid fil- 
tration (to = 0) has a sense only if /ii > 0. This criterion automatically 
implies the equality p2 = oo and that the unique limiting regime for the 
liquid in pores is a rest state. For the case when the crack space is connected 
and pi < oo the author, using the method of reiterated homogenization sug- 
gested by G. Allaire and M. Briane [2], has shown that the limiting velocity 
of the liquid in cracks and the limiting liquid pressure satisfy the usual Darcy 
equations of filtration. For disconnected crack space (isolated cracks), or for 
the case pi = oo the unique limiting regime is a rest state. 

In the present publication we deal with the liquid filtration (tq = 0) and 
the same liquid domain as in Ref . [H] , composed by a periodic system of pores 
with dimensionless size S and a periodic system of cracks with dimensionless 
size e, where S = , r > 1. 

We define the liquid domain Q'y, which is a subdomain of the unit cube 
Q. Let Q = Zf U Zs U 7c, where Zf and Zs are open sets, the common 
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boundary 7c = dZf fl dZg is a Lipschitz continuous surface, and a periodic 
repetition in of the domain Zg is a connected domain with a Lipschitz 
continuous boundary. The elementary cell Zj models a crack space fi^ : the 
domain is an intersection of the cube fl with a periodic repetition in 
of the elementary cell eZf. In the same way we define the pore space flp'. 
Q = Yf UYgU 7p, 7c is a Lipschitz continuous surface, a periodic repetition 
in of the domain Yg is a connected domain with a Lipschitz continuous 
boundary, and flp is an intersection of fl\fll with a periodic repetition in 
of the elementary cell 6Yf. Finally, we put fl'j = fl^ U fl% fl^ = fl\f^ 
is a solid skeleton, and = 9^2^ fl is a "solid skeleton-liquid domain" 
interface. 

Following R. Burridge & J. Keller [8] and E. Sanchez - Palencia [19] we 
describe the joint motion of the mixture of solid and liquid components on 
the microscopic level by well - known system, consisting of the Stokes and 
Lame's equations, coupled with corresponding boundary conditions on the 
common boundary "solid skeleton-liquid domain". For filtration processes 
(tq = 0) we may neglect the inertial terms and consider stationary equations. 
That is, the motion of the incompressible liquid in the liquid domain fl'y is 
governed by the stationary Stokes system 

a^A^- Vg/ + p/F = 0, V-Wf = 0, (0.5) 

for dimensionless microscopic displacements Wf and pressure qf, and the 
motion of the incompressible solid skeleton fl^ is governed by the stationary 
Lame's system 

axAwg - Vqs + PsF = 0, V ■ = 0, (0.6) 

for dimensionless microscopic displacements Wg and pressure g^- On the com- 
mon boundary F^ "solid skeleton-liquid domain" the displacement vectors 
and pressures satisfy the usual continuity conditions 

Wf = Ws, (0.7) 

and the momentum conservation law in the form 

duo 

(a^D(^) - qfl) ■ n = {axH-^s) - qsl) ■ n, (0.8) 

where n{xo) is the unit normal to the boundary at the point Xq G F^. 

In (10. 5p - (10. Sp D(ii) is a symmetric part of the gradient Vu, I is a unit 
tensor, 

2A 

Lgpo 

ps is the mean dimensionless density of the solid phase correlated with the 
mean density of water po and A is the elastic Lame's constant. 
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The problem is endowed with the homogeneous initial and boundary con- 
ditions 

w{x,o) = o, X En = n'jur uni, (0.9) 

w{x, t) = 0, xe S = dQ, t>0, (0.10) 

where w = Wf in Q^j- and w = Wg in Q^. 

Note, that the assumption about incompressibility of the liquid is quite 
natural. It is well - known that the measure of incompressibility is a speed 
of sound of compressible waves. For filtration processes we assume that this 
value is equal to infinity. But the speed of a sound in a solid skeleton in 
two or three times is more than speed of a sound in a liquid. Therefore, we 
may assume that for filtration of incompressible liquid the solid skeleton is 
an incompressible elastic body. 

The case r = 1 corresponds to already studied situation of a simple pore 
space, and the case r > 1 corresponds to a real double-porosity geometry. In 
what follows, we suppose that 

Ho = and < Aq < oo, (0.11) 

where 

lim.a\(e) = An. 

£\0 

For the simple geometry (r = 1) the homogenization procedure has a sense 
only if /xi > (see Ref. [I3]). Moreover, if /xi = oo (extremely viscous 
liquid), then the unique limiting regime is one velocity continuum, describing 
by anisotropic Stokes system for the common velocity in the solid skeleton 
and in the liquid. This fact (that the velocity in the liquid coincides with 
the velocity in the solid skeleton) is a simple consequence of the Friedrichs- 
Poincare inequality. The same situation is repeated for the case r > 1 of 
more complicate geometry. We show that, as before, the homogenization 
procedure has a sense if and only if ni > 0. But this criterion automatically 
implies the equality /i2 = oo. Therefore, due to the same Friedrichs-Poincare 
inequality the limiting velocity of the liquid in pores is proportional to the 
limiting velocity of the solid skeleton. If the crack space is connected and 
/xi < oo, then using the method of reiterated homogenization, suggested by 
G. Allaire and M. Briane [2] we prove that the limiting displacements u of 
the solid skeleton and the limiting liquid pressure g/ satisfy some anisotropic 
Lame's equation 

Ao V ■ (A(") : B{u)) - —Vqf = pF, (0.12) 

m 

coupled with Darcy law for the liquid velocity in cracks 

V, = m,v, + — B(^) [pfF - - Vg/) , (0.13) 
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and common continuity equation: 

V ■ {v,+ {l-m,)vs) =0, (0.14) 

where Vs = du/dt is a velocity of the solid component. 

For the case fii = oo, or for disconnected crack space = ttLcVs and the 
limiting displacements of the solid skeleton and the limiting liquid pressure 
satisfy the usual Stokes system 

Ao V ■ (A(^) : 3(u)) - —Vqf = pF, V ■ u = 0. (0.15) 

m 

Here symmetric and strictly positively definite fourth-rank constant tensor 
A^'') depends only on the geometry of the solid cells Yg and and does not 
depend on criteria Xq and /ii, strictly positively definite constant matrix B^*^) 
depends only on the geometry of the liquid cell Zf and does not depend on 
criteria Aq and fii, p = m pf + {1 — m) ps, m = fy fzX dydz is the porosity 
of the liquid domain, and rric = Xc dz is the porosity of the crack space. 

The system f l0.12p - (10.141) is well - known as Biot's system of poroelasisity 
(Ref. [6]), or Terzaghi system of filtration (Ref. [21]). We call it as Biot - 
Terzaghi system of liquid filtration in poroelastic media. 

Finally, for /ii < oo we consider the family {v^"' '"'^"5 ^f°} solutions 
to the problem (I0.12p - (I0.14p and show that these solutions converge as 
Ao C)0 to the solution of the problem 

V, = — B(^) {pfF - -Vg/) , V-v, = 0, (0.16) 

which is usual Darcy system of filtration and, on the other hand, is a phys- 
ically correct double porosity model for filtration of an incompressible liquid 
in an absolutely rigid body. 

§1. Main results 

To define the generalized solution to the problem (10. 5p - (lO.lOp we char- 
acterize liquid and solid domains using indicator functions in Vt. Let ri[x) 
be the indicator function of the domain VL in M^, that is ri{x) = 1 if a; e 
VL and rjix) = if a; G ]R^\r2. Let also Xpiv) be the 1-periodic exten- 
sion of the indicator function of the domain Yf in Y and Xc{z) be the 1- 
periodic extension of the indicator function of the domain Zf in Z. Then 
Xl{x) = r]{x)xc{x/e) stands for the indicator function of the domain Q^, 
Xp{x) = 77(33) (1 — Xc{x/e))xp{x/S) stands for the indicator function of the 
domain and x^(a3) = xli^) + Xp{^) stands for the indicator function of 
the liquid domain Q"^ . 

We say, that functions {lu^, g^}, where 

= w)x" + <(i - x'), = q}x" + - xl, 
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such that 

w^eL^{iO,Ty,W^in)), ^eL'{iO,Ty,W,'in^j)), q^eL^Gr)) 

is a generahzed solution to the problem (10. 5p - (10.101) . if they satisfy normal- 
ization condition 

/ g^(a;, t)dx = 
Jn 

almost everywhere in (0,T), continuity equation 

\/-w = (1.1) 

in a usual sense almost everywhere in Gt = x (0, T), initial condition (10. 9p . 
and integral identity 

(1.2) 

o 

for any vector-functions if G ^^((0, T)]Wl {Vt)). In (OD 

p' = PfX' + Ps{l - x')- 

The homogeneous boundary condition (lO.lOp is already included into corre- 
sponding functional space. Functions OF /dt and d'^F/dt^ are supposed to 
be - integrable: 

ri = / 1^— I dxdt < oo, r2= I / 1 ^ 1 dxdt < oo. 
Jo Jn Jo Jn 

In the same standard way, as in Ref. pjj, one can show that for any e > 
there exists a unique generalized solution to the problem (10. 50 - (lO.lOp . To 
formulate basic a'priori estimates we need to extend the function w'^ from 
VLI to VL^. To do that we use well-known results (see C. Conca[l0] and E. 
Acerbi at al^) in the following form: for any e > there exists an extension 
w= G L°°((0,T); W^iyt)) such that = in Vt^ and 

!\u'\^dx<Gf \w'\^dx, f p{u')\^dx <G f p{w')\^dx, (1.3) 
Jn Jn% Jn Jn^ 

where G is independent of e and t. 
Holds true 

Lemma 1.1. Let fii > and r > 1. Then there exists sufficiently small 
Eq > 0, such that for any < e < eq and for any < t < T 

\w'{x,t)\'^dx + a^ \I]){w'{x,t))\'^dx + ax \]S){w'{x,t))\'^dx < GF^, 
Jn^j Jn% 

(1.4) 
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\v%x,t)\'^dx + af, \B{v'{x,t))\'^dx + ax \B{v'{x,t))\'^dx < CF2, 



(1.5) 

[ \q%x,t)\^dx= [ {\q}ix,t)\^ + \ql{x,t)\^)dx<C{F, + F2)=CF, (1.6) 
Jn Jn 

^[ \{w' -u')(x,t)\^dx + ^ [ \{w' -u'){x,t)\^dx<CF, (1.7) 

^ ii('''-¥)(--')i^''^^ifii('''-¥)(-')i^''^2 

where = dw'^ /dt and C is independent of e and t. 

Theorem 1. Under conditions ( fO. iij) and conditions of Lemma 2.1 there 
exist functions G L°°((0,T); iy2^(i7)), such that = in a subse- 
quence of small parameters {e > 0}, and functions Vp G L°°((0, T); L^(f2)) - 
the limiting velocity of the liquid in pores, Vc G L°°((0, T); L^(fi)) - the lim- 

o 

iting velocity of the liquid in cracks, u G ((0, T); (fi)) - the limiting 
displacements of the solid skeleton, and qj G L°°((0, T); L^(n)) - the limit- 
ing pressure in the liquid, such that the sequences {Xpd'W^ / dt} , {x'^^^dw'^ / dt} , 
and {gj} converge as e \ weakly in L^((0,T); L^(f2)) to the functions Vp, 
Vc, and qj, respectively. At the same time the sequence {u^} converges as 

o 

e \ weakly in L^((0,T); {VL)) to the function u. 

(I) ///ii = 00, or the crack space is disconnected (isolated cracks), then 

, du du .du du 

Vr, = [i — rrir m„ — — , v.. = rrir — — , v = v^. + v„ + (1 — m) — — = — — , 

and functions u and qj satisfy in Gt the anisotropic Stokes system 

Ao V ■ (A(^) : D(ti)) - — V gf = pF, V ■ u = 0, (1.9) 

m 

with homogeneous initial and boundary conditions 

g/(a;,0) = 0, a; G fi, w(a;,t) = 0, xeS,t>0. (1.10) 

where fourth-rank constant tensor A^^^ is defined below by formula ^3. 30) . 
p = m p f -\- {1 — m) ps , m = jy Jz xdydz - the porosity of the liquid domain, 
nip = Jy Xpdy ^ the porosity of the pore space, and rric = Xcdz - the 
porosity of the crack space. The tensor A^^^ is symmetric, strictly positively 
definite, and depends only on the geometry of the solid cells Yg and Zg. 

(II) If pi < 00, and the crack space is connected, then 

du du 

v„ = (1 — rrir) m„ — — , v = + (1 — rrif)——. 



functions u, v^. and qf satisfy in Gt equations U.9\) . Darcy law in the form 

du 1 1 
v, = m, — + — B(^) {pfF - -Vqj) , x e (1.11) 

initial and boundary conditions U.1(J\) . and boundary condition 

v-n = 0, xeS, (1.12) 

where n is a unit normal vector to the boundary S at x & S . In U.ll\) the 
strictly positively definite constant matrix M^^\ is defined below by formula 
Ii3.18\) and depends only on the geometry of the liquid cell Zj. 

Remark 1.1. Without loss of a generality we may assume that 

/ qj {x, t) dx = 0. 
Jq 

Theorem 2. Under conditions of Theorem 2.1 let yUi < oo and u^''^"\ v^c^^ 
and q^^°^ be a solution to the problem ( (i.gj) Then there exists a 

subsequence of parameters {Xq}, such that the sequence {u^'*'°^} converges 

o 

as Xq y oo strongly in ((0, T); (Q)) to zero, and sequences 

and {qf'^} converge as Xq y oo weakly in L'^{Gt) to functions Vc, and qf 
respectively, which are a solution to the problem 

1 . _ 1 



— b(")| 

jj,i m 



v^= —n--'>[pfF Vg/), xeVL, (1.13) 



V ■ = 0, a; G v^ - n = {), X e S. (1.14) 

§2. Proof of Lemma 1.1 

To prove (11. 4p we choose as a test function in (11.21) the function h{T)dw'^ / dr^x, r), 
where /i(r) = 1, r e (0,t) and /i(r) = 0, r G [t,T): 



f I X^|ro(^(a;,r))pcixcir+^aA j {I - x')H^'{x,t))\''dx 



F ■ ^^^dxdr. 
'0 Jn^ 9t 

Passing the time derivative from dw^ /dt to F in the right-hand side integral, 
applying after that to this integral Holder inequality and the evident estimate 

x'p{w'{x,t))\''dx<C / / x1©(^(a^,r))|2da;cir, 
n Jo Jn c't 
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we arrive at 



J{t) = a^ f x'\Hw%x,t))\^dx + ax f {1 - xn\^iw%x,t))\^dx < (2.1) 

CFi+ f f \w'{x,T)fdxdr. 
Jo Jq 

o 

Next we put Wq = w'^ — . By construction Wq (^/)- To estimate the 
integral 

If= \wlfdx 



we divide it by two parts: 

Pf = ll + Il 4=1 \w^o\'dx, n= I \wl\'dx. 

Let G^^\ where k = {ki,k2,k^) G Z^, be the intersection of fi^ with a set 
{x : x = e{y + k), y e Y}. Then = UkeifiOf^ and 

In each integral 1^ we change variable by x — Sy, then apply the Priedrichs- 
Poincare inequality and finally return to original variables: 

f \w^,\'dx^5' f \wl\'dy< 

5^C^''^ [ \By{wo')\''dy ^ S^C^''') [ \B^{wl)\^dx. 

Here WQ^{y, t) = Wq(x, t), Y^'^^ C F is an appropriate translation to origin of 
the set {l/6)Gi!'\ and C^^) is a constant in the Friedrichs-Poincare inequality 
for the domain y^*"). To estimate these constants uniformly with respect to 
5 (or e) let us clarify the structure of the domain Y^'^\ If the closure of 
has no intersection with the boundary between pore and crack spaces, 
then Y^'"^ = Yj and C*-'^-' coincides with a fixed constant C . Otherwise, 
y*^'^-' is one of two domains, obtained after splitting Y} by some smooth 
surface, asymptotically closed to the plane as £ \ 0. Due to supposition 
on the structure of the solid part Yf, constants C^^^ uniformly bounded for 
all possible planes, splitting Yf. Therefore, supC*^*'^ < C (for simplicity we 
denote all constants independent of e as G) and 

li<5'^Gy^ f \I])^{w'o)\'^dx<6^G [ \B^{w'o)\^dx. (2.2) 
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To explain ideas we consider the easiest geometry, when the hquid part Yf 
is "surrounded" by the sohd part Y^. That is, for each facet S C dY of Y 
the hquid part S fl dYf is completely surrounded by the solid part S fl dY^. 
Due to construction [w^ = in Yg) the constant in the Friedrichs-Poincare 
inequality for Y^''^ depends only on the ratio a = Vj/Vg between the volume 
Vf of the liquid part Y^''^ fl Yj of Y^''^ and the volume K of the solid part 
fl Yg of yC"); C'-'"-' < Ca. It is easy to see, that for chosen geometry 
of Yf and for any type of splitting of Y by planes, this ratio a is uniformly 
bounded. 

In the same way we show that 

I!<e'C [ \B,{wl)\'dx. (2.3) 

In fact, as before we again divide the integral into the sum of integrals 
over domains G^J^^ and make change of variables: 

X = ez, wl{x,t) = WQ{z,t), / {w^l'^dx = / \wl\'^dz. 

For integrals over domains Ci''^ we use the Friedrichs-Poincare inequality, 
based on the fact that the function Wq vanishes on the some periodic (with 
period 6/e) part of the boundary dCi''^ with strictly positive measure, which 
bounded from below independently of e. 
Thus, 

If < 0(5^ + e^) / \B{wl)\^dx < C{— + —)a^ \B{w')\^dx+ 
C{5^ + e^) I p{u')\^dx<C J{t), 



J{t) = a^l p{w^)\^dx + ax I p{w')\^dx, 



and 



/ \w'\'^dx< I \wl\'^dx+ I \u'\'^dx <C{j{t) + ! \w'\^dx). 
Jn'^j Jn^j Jc^ Jni 

To estimate the integral 



\w'^\^dx 

we use the Friedrichs-Poincare inequality, estimate (11. 3p and supposition Aq > 
0: 

I!< [ \u'\^dx<C [ \I]){u')\^dx <Cax [ \B{w')\^dx <CJ{t). 

12 



Gathering all together one has 

/ \w'\^dx < CJ{t). 
Jn 

Estimate (11. 4p follows now from (12.11) and Gronwall's inequality. The 
same estimate (11.41) together with (12.21) and (12.31) result (II. 7p . 

To prove estimates (II. 5p and (11.81) we just repeat all over again for the 
"time derivative" of identity (II. 2p and d'^w^ /dt^. 

Estimate (11.61) is a simple consequence of (ll.4p and (11.51) (see, for example, 
Ref. [13]). 

§3. Proof of Theorem 1 

3.1. Weak and tree-scale limits of sequences of displacements, 
velocities and pressure 

First, we define the velocity of the liquid in pores as = Xpdw'^ /dt, 
the velocity of the liquid in cracks as = xl dw"^ /dt and the velocity of the 
solid skeleton as vl = du^ /dt. By definition 

v^- = vl + vl + {l-x')vl. (3.1) 

On the strength of Lemma 1, the sequences {^'/j, {'V^}^ {'^c}) 

{u^}, {fs}) {Vn^} are bounded in L'^iVtj')- Hence there exists a sub- 
sequence of small parameters {e > 0} and functions qf, Qs, v, Vp, Vc, 

Vs G L^{Gt) and u G L°°((0,T); (fi)) such that 

Qf Qf, qI Qs, ^ V, 
vl Vg, ^ u, Vn^ ^ Vn 

weakly in L'^{Qt) as e \ 0. 
Note also that 

X'a^l^iv') ^ (3.3) 

strongly in L'^{Qt) as e \ 0. 

Next we apply the method of reiterated homogenization (see G. Al- 
laire and M. Briane[2]): there exist functions Q f{x,t,y, z), Qs{x,t,y, z), 
V{x, t, y, z), Vc{x, t, y, z), Uc{x, t, z), and Up{x, t, y, z) that are one-periodic 
in y and z and satisfy the condition that the sequences {o'/}, {qI}) {'*^^}) {'^c}> 
and {Vu^} tree-scale converge (up to some subsequences) to Q f{x,t,y, z), 
Qsix, t, y, z), V{x, t, y, z), V^x, t, y, z), and Vu+VzUdx, t, z)+VyUp{x, t, y, z), 
respectively. The sequence {n^} three-scale converges to the function u{x, t). 

Relabelling if necessary, we assume that the sequences themselves con- 
verge. 



V, 



(3.2) 
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Remind, that three-scale convergence of the sequence {vr^} to the function 
n(a?, t, y, z) means the convergence of integrals 

/ / tt"^ {x,t) ip(^x,t, —, —)dxdt ^ / / / / ^{x,t,y, z) {p{x,t,y, z)dzdydxdt, 
Jo Jn ^ ^ Jo Jn Jy Jz 

for any smooth 1-periodic in y and z function (p{x, t, y, z). By definition the 
function 

n{x,t) = {{Il)y)z, 

where 



(n)y = / Udy, {U)z= I lidz, 
Jy J z 



is a weak hmit in LP'{Gt) of the sequence {tt^}. 
3.2. Macro — and microscopic equations 

We start the proof of the theorem from the macro - and microscopic 
equations related to the liquid motion and to the continuity equation. 

Lemma 3.1. For almost all {x,t) G Gt, y & Y and z & Z, the weak and 
three-scale limits of the sequences {qf}, {ql}, {v"^}, {v^}, {Vp}, and {u'^} 
satisfy the relations 

Qf = —ifix,t)x{y,z), Qs = Qs{i-x), x = Xc{z) + {i ~xc{z))xpiy), 

(3.4) 

Vp = {l-mc)mpVs, V = Vc + {I - mc)Vs, (3.5) 

V-^; = 0, (l-x)(V-w + V. -C/e + Vj, -C/p) =0, (3.6) 

(1 - m)V ■ n + ((1 - x)V. ■ U,)z + (((1 - X)V, ■ UMz = 0, (3.7) 

where m = {{x)y)z - the porosity of the liquid domain, nip = {xp)y - the 
porosity of the pore space, and = {Xc)z ^ the porosity of the crack space. 

Proof. By definition of qj, and and properties of three-scale convergence 
one has equalities Qj = xQfi Qs = (1 — x)Qs- Choosing in (11.21) test 
function in the form = 6 h{t) xj)^ = 5 h{t) i/'(a;, xje^xjb\ where i/?'^ is finite 
in and passing to the limit as e \ we arrive at 

X(y,^) Vj/Q/ = 0, or Qf = x{y,z)Qf{x,t,z). 

Now we repeat all over again with ip = eh{t)'il)^ = e h{t) il){x,x/e), where 
i/?"^ is finite in f2p and get 

x{y,z)V,Qf = 0, or Qf = x{y,z)Qf{x,t), 

which results (13. 4p . 
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(13. 5p is a simple consequence of (13. ip . (11.80 and properties of three-scale 
convergence. 

The first continuity equation in (13. 6p follows from the continuity equation 
(11.11) in the form 

/ v'-V'^dx = 0, (3.8) 
Jq 

which holds true for any smooth functions ip, after passing there to the limit 
as e \ 0. 

Three-scale limit in continuity equation (11.10 in the form 

(1 - x^)V ■ < = 

results the second continuity equation in (13.61) . Finally, (13.71) is just an aver- 
age of the first equation in (13.61) . □ 

Remark 3.2. The first continuity equation in 6]) is understood in the 
sense of distributions as integral identity 

V ■ Vipdx = 0, 

n 

which holds true for any smooth functions ip. 

Lemma 3.2. Let V = {Vc)y- If f^i = oo, then 

V = Vc = Vs{x,t,)xc{z), Vc = mcVs. (3.9) 

If < oo, then for almost every {x,t) G Gt the function V is a 1-periodic 
in z solution to the Stokes system 

-/iiA.y = -V,n--Vg/ + p/F, (3.10) 

m 

V,-V = 0, (3.11) 

in the domain Zj, such that 

V{x,t,z) = v,{x,t), ze^c- (3.12) 

Proof. First of all we derive the continuity equation (13.111) . To do that we 
put ip = eipo{x, x/e) in the integral identity (13. 8p . pass to the limit as e \ 0, 
and get identity 

V ■ V zi>o{x, z)dxdz = 0, 




which is obviously equivalent to (13.111) . 

If /ii = oo, then (13.91) follows from estimate (II. 8p . Let now /ii < oo. 
If we choose in the integral identity (II. 2p a test function in the form 
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= hQ{t)hi{x)il){x/e), where supp /ii C Q, suppt/j(z) C Zj, ■ V' = 0, 
and pass to the hmit as e \ 0, we arrive at 

/ / {h^l^V ■iV,■B,i^i^)) + -qf{\/h■^P)+PJiF■^|^)h^)dxdz = 
Jn JZf 

The desired equation fl3.1Up follows from the last identity, if we pass deriva- 
tives from the test function to V and take into account (13. 111) . The term 
V^n appears due to condition Vz ■ i/' = 0. 

Finally, the boundary condition (I3.12p follows from the representation 

{V)y = V+{l-Xc{z))v,{x,t), 

and inclusion {V)y G W2{Z) for almost every {x,t) G VLt (see Ref. [I3]). □ 

Now we derive macro - and microscopic equations for the solid motion. 
Let ^ ^ 

Qf = ^-(if, Qs = {T-Qs-qf){i-x), qs = {{Qs)zJys- 

Then 

1 1 _ ~ _ _ 

+ Qs) = -^{{Qf + Qs)zs)ys = {{<lf + Qs)z)y = qf + qs 

Ao ^0 

Lemma 3.3. Functions u, U c, Up, qf, and satisfy in Gt the macroscopic 
equation 

V.-((l-m)©H + (l-mp)(©,(C/,))z. + ((0,(t/p))zJn-gl) =F, (3.13) 

where p = m pj + {1 ~ m) p^, q = qf + qs, F = {p/\o)F. 

To prove this lemma we put in (11.21) cp = hQ{t)hi{x), where h is finite in 
Q, and pass to the limit as e \ 0, taking into account (13. 3p . 

Lemma 3.4. Functions u, Uc, Up, and Qs satisfy in Zg and almost every- 
where in Gt the microscopic equation 

V.-((l-Xc)((l-mp)(D(it)+D,(t/,)) + (D,(t/p)-4l)y.)) =0. (3.14) 

To prove lemma we put in (II. 2p (p = ehQ{t)hi{x)(pQ{x/ e) , where hi is 
finite in Q, pass to the limit as e \ 0, and use the equality (1 — x) = 
(1-Xp)(l-Xc). 

Lemma 3.5. Functions u, Uc, Up, and Qs satisfy in Yg and almost every- 
where in Gt x Zg the microscopic equation 

v^- ((i-xp)(m)(^) + ©.(c/c) + %(t/p)-gj)) =0. (3.15) 
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To prove lemma we put in (11. 2p = 5hQ{t)hi{x) ipQ^xje) ifii{x/S), where 
hi is finite in Q, and pass to the hmit as e \ 0. 
3.3. Homogenized equations 

The derivation of homogenized equations is quite standard (see Ref. [13]). 
For the liquid motion we solve the microscopic system (13. 9p - (I3.12p . find V 
as an operator on Vg/ and du/dt, and then use the relation Vc = {V)zf 
Namely, holds true 

Lemma 3.6. Let /ii < oo. Then functions v^, Vg, v = Vc + {I — rric) v^, and 
qj satisfy in the domain Q the usual Darcy system of filtration 

1 1 

= m^v, + — B(')(p/F Vg/), x e Q, (3.16) 

V-i; = 0, xeVL, i;.n = 0, x e S, (3.17) 

where n is a unit normal vector to the boundary S at x E S. 

If the crack space is connected, then the strictly positively definite constant 
matrix B*^'^-' , is defined by formula 

1 3 

B(^) = -V(V^)z,®e,. (3.18) 

In \3.18\) functions V^{z), i = 1,2,3, are solutions to the periodic boundary 
- value problems 

-A,V' + VW = ei, VyV' = 0, zeZf, 
V' = 0, ze 7c, 

where e^, i = 1,2,3, are the standard Cartesian basis vectors and for any 
vectors a, b, and c the matrix a®h is defined as (a® 6) ■ c = a(6 • c). 

If the crack space is disconnected (isolated cracks), then the unique solu- 
tion to the problem [3lM is V' = 0, i = 1, 2, 3, B^'^) = 0, and 

Vc = m^Vs. 

The same procedure is applied for the solid motion. First, we solve the 
microscopic equation (13.151) coupled with the second equation in (13. 6p . find 
Up as an operator on ©^(t/c) and ©(u), and substitute the result into equa- 
tion (I3.14p . Next, we solve the obtained microscopic equation and find Uc 
as an operator on B>{u). Finally, we substitute expressions Up and Uc as 
operators on D(u) into macroscopic equation ( 13.13^ and arrive at desired 
homogenized equation for the function u. 
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Lemma 3.7. For almost every {x,t) G Gt functions u and Uc satisfy in Zg 
the microscopic equation 

V, ■ ((1 - Xc) A(^^ : (OM + e,(C/e))) = 0, (3.20) 

where fourth-rank constant tensor is defined below by formula liS. 23\) . 
Proof. Let 

1 du ' du ' 

As usual, equation fl3.20p follows from the microscopic equations fl3.14p . after 
we insert in the expression 

(©,(t/p))y, - {Qs)ysI = C(^) : (D(n) +D,(C/,)). 

To find it we look for the solution Up to the system of microscopic equations 
fl3.15p and (13.61) in the form 

i,j=l i,j=l 

and arrive at the following periodic boundary - value problems in Yg-. 
Vy ((1 - Xp)((e,(L/;^) + J^) - Qp)) =0, yeY, 

V,- ((i-Xp)(ro.(c/;)-Poi)) =0, yeY, 
V,-t/J + i = o, (t/J)y, = 0, i/en. 

In ((MH) 

jj- = l(p.- + F*) = l(e, ® + e, ® e,). 

Problems (13.211) and (I3.22p are understood in the sense of distributions. 
For example, first equation in (13.2 ip is equivalent to the integral identity 

^(1 - Xp)((D.(C/?) + r^) - Qp) ■■ ^yiv)dy = 
for any smooth and periodic in y function ^p{y). 



(3.21) 



(3.22) 
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The solvability of the problem (13.211) directly follows from the a'priory 
estimate 

/ \VU'^\'dy<C, 

JYs 

and the latter one is a consequence of the energy identity 

I (D,(t/y) : D,(f/;^) + r : ^y{U^^))dy = 0. 

To solve the problem fl3.22p we first find a 1 - periodic function Vq G 
W^{Ys) such that 

Vj,-Vo + l = 0, yeYs. 

There are a lot of ways to construct such a function. In Ref. p!2|, for example, 
one may find non-periodic case. The periodic case is quite similar. 

After that, the solvability of the problem (I3.22p follows from the energy 
equality 

^ (D,(t/;) : (D,(C/;) - By{Vo)))dy = 0, 

which is a result of a substitution into the corresponding to the first equation 
in (I3.22P integral identity the test function {Uq — Vq). 
Thus, 

3 

3 

E ((o,(c/;^'))n - {Q;')Yj^r) ■. (dh + d,(c7,))+ 

((D,(C/;))n ® I - (g;)y. I ® l) : + D,(C/,)) = 

(CS^) + C^^) + + ) : (D(n) + D.(C/c)) = C(P) : (©(-a) + D.(C/j) , 

where B(S)C is a fourth-rank tensor such that its convolution with any matrix 
A is defined by the formula 

(B® C) : A = 1(C : A), 
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and 

3 

= (1 - rrip) r ® + C(^') = (1 - mp)J + C^^), (3.23) 

where 

3 

3 

□ 

Lemma 3.8. Tensors A^^^ and C-^-* are symmetric and the tensor A^'^^ is 
strictly positively definite, that is for any arbitrary symmetric matrices ( = 
(Cij) and rj = {r],j) 

(AW ■.C):V= (A(^) : r/) : C, and (A^^) ■.():(> f3{C : (), 

where positive constant f3 is independent of (. 

Proof. To prove lemma we need some properties of the tensor A*^'^^ which 
follow from equalities 

-{Ql)y^ = {3y{Ul):3y{Ul))Y., (3.24) 
(D,(L7;^-) : D,(C/;))y, = 0, (3.25) 
(Q;^-)y, = -(D,(C/°):FV., (3.26) 

{ByiUin ■■ %(C/?))n + ■■ %(t/?))n = 0, (3.27) 

for all i, j, k, I = 1, 2, 3. 

Equation fl3.24p is a corresponding to the first equation in fl3.22p integral 
identity with the test function Up. Equation fl3.25p is the corresponding 
to the first equation in (I3.22p integral identity with the test function Up . 
Equation (13.260 is the corresponding to the first equation in (I3.2ip integral 
identity with the test function Up. Here we additionally took into account 
relations (13.251) . Finally, equations (13.271) is the corresponding to the first 
equation in (I3.2ip integral identity with the test function [7^'. 



20 



Next we put 

3 3 
i,j=l i,j=l 



Then 

C[^^ : C = mY^))y^, : C = (%(l^?))y., 

and Eqs. (Km - flXTTj) take a form 

(Cf^C):r]=(m).(l^?):©.(l^;))y., (3.28) 

{ByiY,) : D,(l-°))y. = 0, (3.29) 

(C?) : C) : = (C^^ : v) C (3.30) 

(C!^) : r/) : C + (©.(n) ^ %(^.))y. = 0. (3.31) 

Therefore, 

(^) : C) : r/ = (1 - m,)C : V + (C^^^ : C) : = (%(l^;))y. : C+ 



(D,(l-°))y^ : + : (D,(rc))y. + (%(^c) ^ %(^;))i'. + (^ - ^pX ■ V- 
Taking into account fl3.29p and (13.311) we finally get 

(A('=) ■.C):v = (l-m,)C : V+mY^) : By{Y'^))yM%iY'r,))Y. : C+ (3.32) 

(D,(y°))y, : V + (%(n) ■ %(^.))n + C : mY,))YA 
V : (©,(l^c))y. = ((E'.(l"c + + C) : {^yiY, + r^) + r/))y,. 
Eqs. (13.321) and (I3.23P show that tensors A'^'^^ and C^^^ are symmetric: 

In particular, 

(A(^) : C) : C = {{^y{Yc + 1^?) + C) : {%{Yc + 1^?) + C))y. > 0, 

and A'^'^^ is strictly positively definite. In fact, if (A*^^) : C") : C° = for some 
C°, such that (0 : (0 = 1, then 

By{Y^o + X^o) + C° = 0. 

The last equahty is possible if and only if the periodic function l^^o + Y^o 
is a linear one. But due to geometry of the solid cell Ys it is possible only if 
Yt^o + Y'^o = const. Therefore C° = 0, which contradict to supposition. □ 
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Lemma 3.9. Functions u and qj satisfy a.e in Gt the homogenized equation 

V, ■ (a(^) : ©(-a) - p) = |-F, (3.33) 

where fourth-rank constant tensor A^*^ is defined below by formula \3. 36\) . 

Proof. Following the standard scheme, we look for the solution to the micro- 
scopic equation fl3.20p in the form 

3 

where functions U]^ satisfy in Z the periodic boundary - value problem 

V. ■ ((1 - Xc) A(^) : (D.(C/:^) + r^)) = 0, {U'^)z. = 0, (3.34) 
which is understood in the sense of distributions. Thus 

3 

{D^{U,))z. = ( Yl {D.{U'^))z. ® r^) : D{u) = C^^) : D{u), 

3 

C^^^ = J2{D.{Ul^))z.^r', (3.35) 

and 

{{{By{U^) - Qj))n)z. = C(^) : ((1 - me)©(n) + (©,(t/e))zj = 
C(P) : ((1 - me)D(n) + C^^^ : D{u)) = C^p^ : (^((1 - m,) J + C^")) : D{u)^ = 

((1 - m,) C^^') + C^P) : C(")) : D{u), 

A(^) = (1 - m) J + (1 - nip) C(^) + (1 - C^^^ + C^^^ : C^^^^ = 

(1 - m) J + ((1 - mp) J + C^P)) : C^^^ + (1 - m,) C^^^ = 

(1 - m) J + A^") : C^") + (1 - m,) C^^^ = 

(1 - m,) ((1 - mp) J + C^P)) + A(") : C^^^ = 

(1 - rUc) A(^) + A(") : C^") = A^^) : ((1 - m,) J + C^^^) , 

where we have used equalities (1 — m) = (1 — mp)(l — rric) and J : A = A : 
J = A for any fourth-rank tensor A. 
Finally 

A(^) = A^'^) : ((l-mJJ + C('=)), (3.36) 
where C^'^^ is defined by fl3.35p . 

□ 
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Lemma 3.10. The tensor A*-*) is symmetric and strictly positively definite. 
Proof. To prove the second statement of the Lemma we use the equahty 

I (A(") : D.(t/^^')) : ©,(C/f + ! [k^"^ : ©.(J'^')) : B,{U^^)dz = 0, 

J Zs J Z s 

(3.37) 

which is just the corresponding to equation (13.341) integral identity with the 
test function U^^ . 
Let 

3 3 
i,j=l ij=l 

Then (13.371) take a form 

((A(^) : D.(Z^)) : D.(Z^))z. + ((A^^^ : P.(Z^)) : C)z. = 0. (3.38) 
Note also, that by definition 

C('=):C=(D,(Zc))z.. (3.39) 
Relations (I3.38P and (I3.39P result 

(A(^) : C) : = (1 - me) (A^") : C) : ^ + ((A^'^ : C^^^) : c) : ^ = 

(1 - m,)(A(^) : C) : ^ + (A^^^ : (D,(Zc))^J : = (1 - m,)(A(^) : C) : ^+ 
((A(^) : D,(Zc)) : D.(Z,))z, + ((A(^) : D,(Z,)) : C)z. + (A(^) : (©.(^c))^.) ■ V = 

((a^: (D,(Zc)+C)) : (e,(Z,) + r^))^., 
which proves the symmetry of A*^*^ In particular, 

(A(^) ■.rj):rj= ((a^'^) : (D,(Z,) + r/)) : (D,(Z,) + r/))^, > /? (r] : r^). 

□ 

§4. Proof of Theorem 2 

First of all we rewrite the continuity equation in (ILQp and Darcy law 
(II. lip in the form 

V ■ vi^"^ • (B(^) Vq\^°'>) = -pfV ■ (B(")F), (4.1) 

The correctness (uniqueness and existence of the solution) of the problem 
(11.91) - (I1.12P follows from the basic a'priori estimate 

/ / \Vvi^''\x,T)\^dxdT+ — [ \Vqf"\x,t)\^dx<C. (4.2) 
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To derive (14 ■2p we just multiply (14.11) by dq^ ° /dt, and the first equation in 

(11.91) by m vi'^°\ sum results, integrate by parts over domain fi. Integral over 
the boundary S = dfl vanishes due to boundary condition (11.121) . Estimate 
(14.21) follows now from Holder, Gronwall and Korn's inequalities. Next we ap- 
ply the standard compactness results to choose the convergent subsequences 
of {vi^°^} and {q^f^°^}, and pass to the limit as Aq /" oo in (II. lip and in the 
integral identity, corresponding to the continuity equation in (II. 9p . Estimate 
(14.21) also guarantees the strong convergence of {vi'^°^} to zero as Aq /" oo. 

Conclusions 

We have shown how the new rigorous homogenization methods can be used 
to clarify the structure of mathematical models for liquid filtration in natural 
reservoirs with very complicate geometry. Obvious advantage of suggested 
models are: 

1) their solid physical and mathematical bases - the models are asymp- 
totically closed to trustable mathematical model on the microscopic level; 

2) their clear physical meaning - the choice of the model depends on 
ratios between physical parameters of a process in consideration; 

3) for most often met situation of disconnected crack space the suggested 
model is so simple as well as usual Darcy system of filtration, but, in contrast 
to the last one, its solutions are more regular, that is very important in 
applications to various nonlinear problems. For example, at the description 
of replacement of oil by water. 
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